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ABSTRACT 


The multilayer theory of anisotropic elasticity and a finite element method were used to analyze the open-mode 
delamination stress concentrations in horseshoe and elliptic laminated composite curved bars. Two types of lamina- 
tions, solid laminations and sandwich laminations, were analyzed. 

It was found that the open-mode delamination stress concentration could be greatly increased in these two types 
of curved bars by decreasing their aspect ratios. The open-mode delamination stress concentration generated in the 
solid laminations was found to be far more severe than that generated in the sandwich laminations. 

The horseshoe curved bar may be used to determine both the open-mode delamination strength of solidly lami- 
nated composites and the open-mode debonding strength of sandwiched laminated composites. However, the elliptic 
curved bar is only good for determining the open-mode delamination strength of solidly laminated composites. 


1 INTRODUCTION 

One of the shortcomings of the laminated composite structures is the delaminations between the composite 
laminas. Once the delaminations take place, the composite structures could lose their structural integrity (that is, 
stiffness, service life, etc.) and could result in ultimate catastrophic delamination failure. The failure could be either 
“open-mode” or “shear-mode” delamination failure depending on the stress field induced in the composite structures. 

In the applications of the composite materials, certain structural components could contain curved beam regions, 
or the entire structural components could be in the form of curved panels (either solid laminations or sandwich 
laminations). When such structural components are subjected to cyclic bendings during service loadings, both radial 
(thickness direction) and shear stresses could be induced in the structures. These two types of stresses are distributed 
nonuniformly and their peak values (stress concentration points) are located at different sites where delaminations 
between composite layers or between face-sheets and sandwich core (debonding) could nucleate. The open-mode 
delamination can occur only when the bending is in the direction of increasing the curvatures of the curved structural 
components. However, the shear-mode delamination is invariant to bending directions. 

In measuring the open-mode delamination strength of the composite structures, usually a flat square specimen 
is pulled in the thickness direction (see fig. 1(a)). In this type of test, debonding between the specimen and the load 
application fixtures (platters) can occur instead of composite delamination. In addition, it is difficult to maintain 
perfectly parallel separation of the two platters during this type of test. An alternative geometry for the open-mode 
delamination test specimen is in the shape of a semicircular curved bar shown in figure 1(b). In this type of test 
specimen, the open-mode delamination can nucleate at the peak tensile radial stress point located at the midspan 
of the curved bar. Through the use of the theory of anisotropic elasticity, it is possible to calculate the open-mode 
delamination strength of the composite without direct measurement at the delamination site. The shortcoming of the 
semicircular curved bar specimen is that under the end forces, the magnitude of the open-mode delamination stress 
is identical to that of the shear-mode delamination stress Oocated at either end of the specimen) (ref. 1). Thus, if the 
composite has weaker delamination strength in shear mode than in open mode, the open-mode delamination cannot 
be achieved in the semicircular curved bar specimen. The way to circumvent this problem is to generate open-mode 
delamination stress concentration by changing the specimen geometry to horseshoe shaped or elliptic shaped, etc. 
This paper presents the studies of the open-mode delamination stress concentrations which could be generated in the 
said two types of test specimens fabricated with solid laminations or sandwich laminations. 
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arbitrary constants associated with nonpure bending 
inner radius of semicircular region of horseshoe curved bar 
outer radius of ith layer of semicircular region of horseshoe curved bar 
mean radius of semicircular region of horseshoe curved bar, j( a + b) 
arbitrary constants associated with pure bending 

outer radius of semicircular region of horseshoe curved bar, or semi-minor axis of 
outer boundary of elliptic curved bar 

semi-major axis of outer boundary of elliptic curved bar 

modulus of elasticity of single ply in fiber direction 

modulus of elasticity of single ply transverse to fiber direction 

modulus of elasticity in r direction 

modulus of elasticity in 0 direction 

quadrilateral membrane element 

moment arm, M = Pe 

shear modulus of single ply 

shear modulus associated with r-6 system 

width of sem icircular curved bar 

index associated with ith layer, i = 1,2, 3,. ..,7V 

open-mode stress concentration factor, [( <r r ) max of any curved bar]/[( a r ) max of 
semicircular curved bar] 

anisotropic parameter, associated with pure bending 

applied end moment, M - Pe 

total number of laminated layers 

applied end force 

radial distance 

radial location of zero <j$ 

structural performance and resizing 

thickness of curved bar 

sandwich core thickness (or depth) 

thickness of sandwich face sheets 

displacement in r-direction 

displacement in 0-direction 

rectangular Cartesian coordinates 

anisotropic parameter, 1—2 v r e) + Jk- associated with nonpure bending 

tangential coordinate 

Poisson ratio of single ply composite 
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Poisson ratios associated with r-6 coordinates 

distance in thickness direction measured from inner boundary of elliptic curved bar 
radial stress 

debonding stress calculated from strength of materials 
open-mode delamination stress (maximum tensile radial stress) 
radial stress induced by the end moments M 
radial stress induced by the end forces P 
tangential stress 
shear stress 

quantity associated with ith layer 
quantity associated with ith layer 

2 COMPOSITE CURVED BARS 

2.1 Horseshoe Curved Bar 

Figure 2 shows the geometry of the horseshoe curved bar which consists of a semicircular region and two straight 
regions at both ends. This type of curved bar specimen could be ideal for testing the open-mode delamination strength 
of the solidly laminated composites or sandwich laminated composites. 

By properly adjusting the end moment arm e (the length of the end straight region), it is possible to increase 
the intensity of the peak radial stress several times that of the peak shear stress. Thus, the open-mode delamination 
failure can be achieved (instead of delamination failure in shear). 

The bending of the horseshoe curved bar under end forces P can be analyzed by the superposition of the following 
two cases shown in figure 3 (for small deformations): 

1. Semicircular curved bar subjected to end forces P (nonpure bending), 

2. semicircular curved bar subjected to end moments M - Pe (pure bending). 

By using the multilayer theory developed by Ko (ref. 2) for the pure and nonpure bendings of a laminated semi- 
circular composite curved bar (solid or sandwich laminations), it is possible to calculate the open-mode delamination 
stress in terms of the applied load (measured) without direct measurement (which is nearly impossible) of the de- 
lamination stress at the delamination nucleation site. In the analysis the thicknesses of the two composite face sheets 
of the sandwich lamination will be set equal (see fig. 2). The open-mode delamination stress will be calculated for 
different aspect ratios of the horseshoe curved bar for studying the effect of the aspect ration on the delamination 
stress concentration. 

2.2 Elliptic Curved Bar 

The second type of curved bar (candidate geometry for the open-mode delamination test specimen) analyzed 
was a semielliptic curved bar shown in figure 4. The outer boundary of this curved bar is semielliptic, but the inner 
boundary is not elliptic due to constant thickness. No attempt was made to develop a multilayer theory for this type 
of curved bar. Instead, a finite element method was used to analyze the open-mode delamination stress concentration 
in this type of curved bar (solid and sandwich laminations) for different curved bar aspect ratios for comparison with 
the case of horseshoe curved bar. 
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3 MULTILAYER THEORY 


The multilayer theory developed by Ko (ref. 2) (using Lekhnitskii’s anisotropic stress functions (ref. 3)) for a 
laminated anisotropic semicircular curved bar under pure or nonpure bending may be used to calculate the open- 
mode delamination stress for the horseshoe curved bar. If the semicircular curved bar is fabricated with N number 
of anisotropic lamina, the stresses and the displacements induced in the ith layer (i = 1 , 2 , 3 , . . . , N) can be written 
as follows (ref. 2). 

3.1 Stresses 


a. For end forces P 


C 4 l) (r,0 ) = 
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b. For end moments M 


a^ir) - 2 B'i+ C'( 1 + kp r ki ~ 1 + D\{ 1 - kp r~ k '~ x (4) 

(r) = 2B[+ C’M 1 + *i) r ki ~ l - D’M 1 - *i) r (5) 
= 0 ( 6 ) 

where Ai, Bi,Di and B[, C[, D'i are two sets of unknown arbitrary constants (associated with the ith layer) which 
are determined from boundary conditions shown in section 3. Notice that the magnitudes of (eq. (1)) and r^g 
(eq. (3)) for the case of end forces P are identical, but are out of phase in the 9 direction by j. 

3.2 Displacements 
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b. For end moments M 
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where the rigid body motion terms are neglected. 

3.3 Determinations of Arbitrary Constants 


k.E^ 


JO ' 

r?(0 


(9) 

( 10 ) 


To determine 3 N unknown constants A,-, B i} A or C< , D\ ( i = 1 , 2 , 3 , . . . , N) , 3 N boundary conditions 
are required. Those boundary conditions are listed in the following. 

a. For end forces P 

The traction-free conditions at the inner boundary (r = a) and at the outer boundary (r = 6) give two conditions: 


v { r l) (a,0) = 0 ( 11 ) 

4 N) (b,0) = 0 (12) 

And the continuities of radial stress cr r , radial displacement u r and tangential displacement uq at each interface 
between layer i and i + 1 (i — 1,2, 3,..., AT — 1) for perfect bonding give 3(N — 1) conditions: 

= a< i+1) (ai,0) (13) 

u< o (ai,0) = u[ i+X \a u e) (14) 

*4 l) (ai,0) = Uq 1+1 \ 0,1,0) (15) 


Notice that in the traction-free conditions at the free surfaces, and in the stress continuity conditions at each 
interface, the shear stress t t q was not included because it has identical r-dependency as <j r (see eqs. (1) and (3)). 
Lastly, the expression of the applied force P in terms of the shear stress t t q at the end (9 = 0) of the curved bar 
gives one condition: 

-P = Y, r T re(E0)dT (16) 

i= 1 ^ a *- 1 

Equations (11) through (16) will give 3 N conditions for the determination of 3 N unknown constants Ai, £?i, D% by 
solving 3N x3N simultaneous equations numerically. 

b. For end moments M 

Similarly the traction-free conditions and the interface stress and displacements continuity conditions give the 
following 3iV — 1 conditions: 


o 

II 

b 

(17) 

o { r N) (b) = 0 

(18) 

cr£ l) (ai) = a^ ,+ 1) (oi) 

(19) 

4 l ' } (a;) = 4 ,+ 1) (a,) 

(20) 

^(a,-) = M0 1+1) (fli) 

(21) 
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The relationship between the applied moment M and the tangential stress <jq gives the N\h condition: 

* r ai 

= (r — ro)ce(r)dr (22) 

i=i ^ a *- 1 

for the determination of 3 N unknown constants C [ , D[ by solving 3 Nx3N simultaneous equations numerically. 
The ro appearing in equation (22) is the unknown radial location where <jq is zero. For pure bending the integration 
of the ro term vanishes. 

3.4 Superposition of Solutions 

The two stress fields calculated for the semicircular curved bar under end forces P (nonpure bending) and end 
moments M (pure bending), may be superimposed to generate a stress field for the horseshoe curved bar (see fig. 3). 
The open-mode delamination stress (<J T )max (located on the z-axis, fig. 2) will then be the maximum value of the 
summation of two radial stresses a^( , f ) due to P and o^( vp) due to M 

(ov)max = of (r D , 0 + a^(r D ) (23) 

where rjj is the radial location of (cr r ) maa; , and is determined from the condition ^ = 0 which gives 

AiPiift - 1 )r§ + BiPiifc + 1 )r~ I f i - D { + (kj - 1 )(C<r$ + D\r~ D k ') = 0 (24) 

assuming ( a T ) max occurs in the ith layer. 

4 FINITE ELEMENT ANALYSIS 

Finite element stress analysis was performed on the multilayered elliptic curved bar (see fig. 4) by using the 
structural performance and resizing (SPAR) finite element computer program (ref. 4). Figure 5 shows the SPAR finite 
element model set up for the elliptic curved bar (solid lamination or sandwich lamination). Because of symmetry 
with respect to the x-axis, only half span of the curved bar was modeled. The elliptic curved bar was modeled with 
22 layers of anisotropic quadrilateral membrane E41 elements in the thickness direction, and 90 E41 elements in the 
circumferential direction. For the case of sandwich lamination, each of the face sheets was modeled with 5 layers 
of E41 elements in the thickness direction, and the sandwich core with 12 layers of E41 elements in the thickness 
direction. The points in the midspan (along x-axis) were allowed to move freely only in the x -direction. This is 
represented in figure 5 by a train of rollers. At the upper end of the model, only one point lying in the middle surface 
was constrained in such a way that it could move freely only in the y-direction. This is shown in figure 5 with only 
one roller. Thus, the end of the model is allowed to rotate freely. 

5 NUMERICAL EXAMPLES 

The type of composite used in the analysis has the following ply material properties 

El = 17.2369 x 10 10 N/m 2 (25 x 10 6 lb/in 2 ) 

E t = 0.8274 x 10 10 N/m 2 ( 1.2 x 10 6 lb/in 2 ) 

Glt = 0.4137 x 10 10 N/m 2 (0 .6 x 10 6 lb/in 2 ) 
vlt ~ 0 -33 
v TL = 0.01584 
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which give the values of anisotropic parameters as (3 = 7 .9272 and k = 4 .5644. 

The honeycomb core used in the sandwich laminates has the following material properties 

E r = 6.8948 x 10 8 N/m 2 (0 .1 x 10 6 lb/in 2 ) 

Eg = 6.8948 x 10 4 N/m 2 ( 10 lb/in 2 ) 

G r g= 1.2411 x 10 8 N/m 2 (0.018 x 10 6 lb/in 2 ) 

Vg r — 6.5 X 10- 5 

The semi-minor axes b of both solid and sandwich laminations for the two curved bar geometries (see figs. 2 and 
4) were set to be b = 5 .08 cm (2 in.), but the semi-major axes c of both types of laminates were allowed to change 
for different aspect ratios. 

All the solid laminations were fabricated with 22 layers of 0° -plies (fibers in the circumferential direction), 
having the total thickness of t = 0 .3175 cm (0.125 in.). All the sandwich laminations have the total thickness of 
t - 1 .1608 cm (0.457 in.), with a sandwich core thickness of t c = 0 .6350 cm (0.25 in.), and a face sheet thickness 
of t f = 0.2629 cm (0.1035 in.). Each of the face sheets contains 18 layers of 0° -plies. This type of sandwich 
lamination was used in the fabrication of X-wing helicopter rotor blades. For a special case of aspect ratio b/c — 1 
(that is, a semicircular sandwich curved bar), a sandwich core of ±15° lamination was also studied. This type of 
core has the anisotropic parameters of (3 = 4 .2498 , k = 2 .4304. 

6 RESULTS 

6.1 Solid Laminations 

Figure 6 shows the distributions of tensile radial stress cr r along the x-axis for different values of aspect ratio b/c 
of the horseshoe curved bar with solid lamination. The open-mode delamination stress ( ov) max (maximum value of 
o>) could be increased many times that of the semicircular case (~ = 1 ) by lowering the aspect ratio b/c of the curved 
bar. For example, the open-mode delamination stress concentration factor K (( <J r )max of any \ value normalized by 
((j r )max of semicircular case (~ = 1)) for ^ = 0.1 reached as high as 14.4431 (for the anisotropic case). Notice that 
as the value of b/c decreases, the location of (<J r )max point migrates toward the inner boundary of the curved bar. 
In figure 6, the stress curves for the isotropic case ((3 = 2 , k = 1 ) were also shown for comparison. The difference 
between the anisotropic and isotropic cases became more conspicuous for lower aspect ratios ^<0.2. The isotropic 
case gave slightly higher value of ( & r )max- Figure 7 shows the plots of cr r along with the z-axis for different values 
of b/c of the elliptic curved bar fabricated with solid lamination. For the elliptic curved bar, the rate of the increase 
of K with the decrease of b/c is faster than that of the horseshoe curved bar (see fig. 8). The reason is that for the 
same aspect ratio b/c , the former has smaller midspan radius of curvature than that of the latter. Notice from figure 8 
that the value of K given by elliptic curved bar at £ = 0 .3 is slightly higher than that given by the horseshoe curved 
bar at £ = 0 .1 (for the anisotropic case). 

6.2 Sandwich Laminations 

Figure 9 shows the distributions of <j r along the x-axis for various aspect ratios of the horseshoe sandwich curved 
bar. Unlike the case of solid lamination (see fig. 6), ( <J r )max occurred in the inner face sheet for the aspect ratios 
- c > 0 .3. For \ < 0 .2 , the location of ( a r )max moved to the interface between the inner face sheet and the sandwich 
core. For testing the debonding strength of the sandwich laminations, the aspect ratio - c of the horseshoe sandwich 
curved bar test specimen should be designed for ^<0.2. 
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Figure 10 presents the distributions of cr r along the rr-axis for different values of b/c of the elliptic sandwich 
curved bar. For this case, ( cr r ) max is located within the inner face sheet for the whole range of aspect ratios. Thus 
the elliptic curved bar specimen is not recommended for testing the debonding strength of the sandwich lamination 
because the interfacial debonding failure between the face sheet and the sandwich core may not occur if the inter- 
laminar delamination strength of the face sheet is weak and the specimen could fail under interlaminar delamination 
of the inner face sheet. Again the stress distributions in the sandwich core for all aspect ratios are nearly linear. 
Figure 11 shows the comparison of the rate of increase of (c y T )max with the decrease in b/c for the horseshoe and 
elliptic sandwich curved bars. The magnitude of ( (J r ) m a X for the sandwich laminations is somewhat lower than that 
of the solid laminations (see fig. 8). At £ = 0 .4 , the elliptic case gave the value of ( a r ) max as being slightly higher 
than that given by the horseshoe case at £ = 0 . 1. 

6.3 Semicircular Sandwiches 

As presented in section 6.2, for honeycomb core sandwich laminations, ( (J r ) max never occurred within the hon- 
eycomb core. The study of the semicircular sandwich curved bar is to show how the core material properties will 
affect the location of (oy) ma:r . Figure 12 shows cr r distribution in the 9 = j plane for three cases of sandwich core 
of the semicircular sandwich curved bar under end forces P (nonpure bending). For the honeycomb core (solid 
curve), the (a r ) max is located in the inner face sheet. As the core material becomes stiffer (±15 lamination, long- 
short broken curve), the ( <j T ) max location migrates into the core. The curve for solid lamination is also shown for 
comparison. Figure 13 shows the case of pure bending of the semicircular sandwich curved bar. Like the nonpure 
bending, the stiffer core of ±15 lamination caused the location of (<J r ) max to move within the core (long-short 
broken curve). For the honeycomb core (solid curve), (cr r ) maa: is always located at the interface between the inner 
face sheet and honeycomb core. The debonding stress cr* calculated from the strength-of-materials equation (see 
fig. 14) 

* = M_ 

<Jr hta 

is also shown for comparison. The value of cr* is slightly lower than ( <j r ) max calculated from the multilayer theory 
(solid curve) and finite element method (solid dots). 

7 CONCLUDING REMARKS 

The multilayer theory of anisotropic elasticity and a finite element method were used to analyze the open-mode 
delamination stress concentrations in horseshoe and elliptic laminated composite curved bars of different aspect 
ratios. Both the solid laminations and sandwich laminations were analyzed. The open-mode delamination stress 
concentration can be greatly increased in these two types of curved bars by decreasing their aspect ratios. The open- 
mode delamination stress concentration generated in the solid lamination was more severe than that generated in the 
sandwich lamination having the same aspect ratio. The elliptic curved bar was slightly more efficient in generating 
the open-mode delamination stress concentration as compared with the horseshoe curved bar. The horseshoe curved 
bar could be used for the determination of both the open-mode delamination strength of laminated composites and 
the debonding strength of the sandwiched composites. The elliptic curved bar is not suitable for determining the 
debonding strength of the sandwich composites. The location of the open-mode delamination stress is a function of 
the curved bar geometry and the material properties. 


Ames Research Center 
Dryden Flight Research Facility 
National Aeronautics and Space Administration 
Edwards , California , September 26, 1989 
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Figure 1. Geometries of composite delamination test specimens. 
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Figure 8. Open-mode delamination stresses in curved bars as functions of aspect ratio 







Figure 12. Distributions of radial stress in 9 = 7t/2 plane of a semicircular sandwich curved bar under end forces. 
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Figure 13. Distributions of radial stress in a semicircular sandwich curved bar under end moments. 
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types of laminations, solid laminations and sandwich laminations, were analyzed. 

It was found that the open-mode delamination stress concentration could be greatly increased in these two 
types of curved bars by decreasing their aspect ratios. The open-mode delamination stress concentration 
generated in the solid laminations was found to be far more severe than that generated in the sandwich 
laminations. 

The horseshoe curved bar may be used to determine both the open-mode delamination strength of solidly 
laminated composites and the open-mode debonding strength of sandwiched laminated composites. How- 
ever, the elliptic curved bar is only good for determining the open-mode delamination strength of solidly 
laminated composites. 
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